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We present a detailed analysis of a quantum memory for photons based on controlled and reversible 
inhomogeneous broadening (CRIB). The explicit solution of the equations of motion is obtained in 
the weak excitation regime, making it possible to gain insight into the dependence of the memory 
efficiency on the optical depth, and on the width and shape of the atomic spectral distributions. We 
also study a simplified memory protocol which does not require any optical control fields. 
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I. INTRODUCTION 

The implementation of quantum memories for photons 
is an important goal in quantum information processing. 
It would allow the realization of on-demand single-photon 
sources based on heralded sources [l|, and provide a ba- 
sic ingredient for quantum repeaters [3| . Several different 
approaches to the realization of such memories have been 
proposed, using both single absorbers in high-finesse cav- 
ities [3] and dense atomic ensembles. The latter proposals 
include the use of off-resonant interactions ^ , of electro- 
magnetically induced transparency (EIT) [6], and of non- 
standard photon echoes On the experimental side, 
storage and retrieval of classical light has been realized 
using EIT ^ and photon echoes [9j. Coherent states of 
light have been stored in an atomic ensemble with higher- 
than-classical fidelity using off- resonant interactions jlfl] . 
Recently the storage and retrieval of single photons have 
been reported using EIT in atomic ensembles [HI, [T^ • 

In the present work we are following the approach that 
originated in Ref. Q , where it was shown that highly ef- 
ficient photon echoes [l^ can be generated in a gas of 
atoms exploiting the fact that the Doppler shift changes 
sign if the propagation direction of the light is reversed. 
A first adaptation of the effect to solid state systems 
was su gge sted in 114] using nuclear magnetic resonance. 
Refs. and proposed an attractive experimen- 

tal realization of the same principle using controlled re- 
versible inhomogeneous broadening (CRIB). These pro- 
posals combine spectral hole burning techniques and the 
use of controllable Stark shifts from electric field gra- 
dients. First proof-of-principle experimental demonstra- 
tions of the CRIB approach have recently been performed 

[Sill- 

Here we perform a detailed theoretical analysis of the 
CRIB quantum memory protocol for finite optical depth 
and general atomic distributions. In section[TT]we give the 
equations of motion for the full system atoms plus light 



and recall the principle of a memory based on CRIB. In 
section IIIII we derive the general solution of the equa- 
tions of motion under the condition of weak excitation. 
We discuss both the complete memory protocol, where 
the output field is emitted in the backward direction, and 
a simplified protocol, which does not use any optical con- 
trol fields, leading to forward emission. In section HVl we 
study the dependence of the memory efficiency on the op- 
tical depth of the medium. We show that it approaches 
one for the complete protocol for sufficiently large optical 
depth, whereas it can be over 50 % for the simplified pro- 
tocol, for which the efficiency is limited by reabsorption. 
In section |V] we show that the width of the initial atomic 
distribution determines the possible storage time. We 
furthermore analyze what is the optimal broadening for 
a given initial distribution and pulse width. Finally we 
investigate the effect of shapes of the atomic distribution 
and light pulse. Section IVTl gives our conclusions. 



II. GENERAL PRINCIPLE 

A light pulse propagates through a medium composed 
of two-level atoms. It resonantly couples the two atomic 
states, the ground state \g) and the excited state |e). We 
look at the behavior of the positive frequency part of 
the slowly time- varying envelope E{z, t) of the light field 
decomposed in forward and backward modes 

E{z, t) = E{{z, <)e*""^/'= + Ebiz, t)e~"^°'/'' (1) 

in a onc-dimensional light propagation model. This one- 
dimensional model is well adapted to the propagation 
either in an optical wave guide or in a bulk medium under 
the condition that A/X£ ^ 1, A being the area of the light 
beam, A the wavelength and i the propagation length. 
The central frequency ujq of the light pulse is detuned 
from the atomic transition frequency LUcg by A cjo — 
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The atomic transitions undergo an inhomogeneous 
broadening that we consider uniform in space such that 
the density of atoms associated to the detuning A is p(A) 
p^ . To describe the properties of the atomic ensemble. 
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we define a mean field per atoms, slowly varying in time 

i,j standing for e or g. In the above sum, the atom in- 
dex n runs over all atoms N{A, z) := p{A)dz5A/L with 
detuning within the interval [A — SA/2, A + 5A/2] and 
position within the interval [z — dz/2, z + Sz/2]. L denotes 
the length of the medium. The Hamiltonian describing 
the system atoms plus light is given in the rotating wave 
approximation by 

H = r dA^ I dz[Aaee(z,t;A) (3) 

J-oo ^ Jo 

-pE{z,t)aeg{z,t;A) + H.c.] 

p being the dipole moment of the transition |e)-|(7). In 
analogy with the light field, the positive frequency part 
of the operator associated to the atomic coherence ago is 
decomposed into two counter-propagating contributions 

age(z, t; A) = m{z, t; A)e''^°'^' +a^{z, t; A)e-^""^/^ (4) 

Under the approximation that most of the population 
stays in the ground state dgg « 1, which is well justi- 
fied for single or few photons light storage, the evolu- 
tion of the atomic coherence is given by the following 
Heisenberg-Langevin equations 

d 

— (Tf (z, t- A) = AcTf (z, t- A) + ipEf{z, t), (5a) 



— (JbC^:, t; A) = -iAcrb(z, t; A) + ipEh{z, t). (5b) 
ot 

Here we neglect the homogeneous decoherence. We will 
however take into account inhomogcneous dephasing in 
the following, cf. section IV Al Using the slowly varying 
approximation, the evolution of the forward and back- 
ward components of the light pulse is given by 

[dt + ""YzJ j ^A G(A)af (z, t- A), 

(6a) 

[j^-c—\E^{z,t)=i(3 J dAG(A)ab(z,i;A) 

m 

where (3 is defined hy (3 :— g^Np with go :— \J bjQ / {2€qV) . 
N := J dAp{A) is the number of atoms in the quantiza- 
tion volume V such that p(A) = iV x G{A) where G(A) 
is the normalized spectral atomic distribution. 
A complete description of the system is given by the set of 
Eqns. (O-®. These equations of motion describe a sit- 
uation in which the presence of the laser field introduces 
changes in the medium by inducing a polarization of the 



atomic ensemble. Reciprocally, the atomic polarization 
can serve as a source for the optical field. Owing to 
the consideration of weak excitations, the coupled equa- 
tions are linear and of first order offering the possibility 
to solve them analytically. Furthermore, the linearity of 
the equations of motion implies the equivalence between 
classical and single-photon dynamics. In the following, 
the quantities Ey, and E{ can be interpreted either as 
classical fields or as single-photon wave-functions. Anal- 
ogously, CTb and (Tf can describe either a classical atomic 
polarization or the wave-function of a single atomic ex- 
citation created by a single photon. 

The principle of a memory based on reversible absorption 
can be deduced from the symmetry analysis of the equa- 
tions of motion. The forward components of the light 
operator and of the atomic coherence evolve following 
the set of equations (|5ap - ((6a|) . If the following transfor- 
mations are performed 

A ^ -A, (7a) 

El, ^ -Eb, (7b) 

the equations for the backward components (j5bp - (j6bp de- 
fine a time reversal evolution of the system compared to 
the equations for the forward components. This analysis 
reveals the essence of a quantum memory based on CRIB 
,7,, .16]: a light pulse propagating in the forward direction 
is completely absorbed by the atomic ensemble. What 
is required to retrieve the pulse as a time reversed copy 
of itself is to reverse the detuning between the spectral 
components of the light pulse and the atomic transition 
frequency. At the same time, one has to apply a phase 
matching operation such that the wave function of the 
single atomic excitation propagates in the backward di- 
rection. 

We now describe how one can construct an optical mem- 
ory according to the principle discussed above. The 
followin g sc enario closely follows the ideas presented in 
Refs. [Elll. 

(i) One first chooses the medium. The relevant atomic 
transition should have a frequency in the optical domain. 
It should also have low decoherence rate compared to the 
total duration of the pulse in order to keep the phase 
and amplitude informations during the optical absorp- 
tion. In Ref. (l6j . several experimental realizations are 
discussed using atomic gases at room temperature or ion- 
doped solid-state materials cooled to cryogenic tempera- 
ture. Motivated by the work done in our group f20l. l2ll| . 
we here consider a medium based on ion-doped solid ma- 
terials. 

(ii) A narrow atomic absorption line is prepared. In 
ion-doped materials, the absorption transition is broad- 
ened inhomogeneously due to the environment in the host 
which is different for each ion. A narrow absorption line 
can be prepared using optical pumping techniques. The 
width of the prepared line is limited ultimately by the 
homogeneous line-width defined from the decoherence 
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FIG. 1: (Color online) Schematic representation of the spec- 
tral atomic distribution. The initial distribution Go(Ao) with 
characteristic bandwidth 70 is represented as black dotted 
line. Three spectral components of the initial distribution 
are considered (blue, black and red vertical lines) and each of 
them is broadened according to a distribution G'(A') (blue, 
black and red dashed lines) with bandwidth 7. The final 
broadened distribution, called G(A) is the convolution of the 
initial distribution Go(Ao) and of the broadened distribution 
G'(A') associated to a single initial absorption line (Ao). The 
pulse shape with bandwidth T is represented as green dashed- 
dotted line. 



rate of the atomic ensemble. In practice however, the 
bandwidth of the absorption line might be limited by 
the bandwidth of the laser used for the optical pumping. 
We will denote the initial detuning by Aq. The selected 
inhomogeneous atomic distribution, called G'o(Ao) with 
bandwidth 70 is represented in Fig. [l]as dotted line, 
(iii) The initial distribution Go(Ao) associated to the 
detunings Aq is then broadened to a larger bandwidth 
using reversible and controlled inhomogeneous broaden- 
ing. Every detuning Aq is broadened to Aq + A', where 
A' is distributed according to G"(A') (dashed line in 
Fig. [1]). In ion-doped materials, a controlled inhomo- 
geneous broadening of the initial absorption line can be 
obtained by applying an electric field gradient that shifts 
the transition frequency of the ions by the Stark effect 
such that the detuning A' is given by the induced Stark 
shift [ll,[i3,[il[20|. The final distribution G(A) (fuU line 
in Fig. [T]), taking into account the broadening of all the 
detunings Aq, is given by the convolution of the initial 
distribution G'o(Ao) and of the distribution G"(A'). By 
reversing the electric field, the Stark shift is reversed and 
consequently, each detuning is reversed such that Eqn. 
([7a|) is replaced by 



Ao -I- A' ^ Ao - A'. 



(8) 



Here we consider a controlled inhomogeneous broadening 
which is independent of the position z. This corresponds 
to a rare-earth doped fiber or to a crystal in which the 
electric field gradient is applied in a direction transverse 
to the propagating pulse. 

(iv) The equations of motion associated to the forward 
components are connected to the ones associated to the 
backward components by applying a position-dependent 



phase e^*""^/"^ to the atomic ensemble when the light 
excitation is reduced to zero. This transforms the 
forward component df into the backward component 
ab according to the Eqn. (g]). The phase shift e^*""^/'' 
can be realized by using counter-propagating 7r-pulses 
to transfer the atomic coherence back and forth to an 
auxiliary ground state p^ . [iGj . This procedure can 
also offer longer storage duration by choosing a ground 
state with long coherence time [131 . If the phase shift 
is not applied, the evolution of the system is given by 
the forward set of equations. The pulse reemitted by 
the medium propagates in the forward direction and 
is partially reabsorbed, limiting the memory efficiency 
On the other hand, this simplified protocol does not 
require the use of additional control laser fields. This 
provides an advantage for few photons quantum memory 
by avoiding the light noise generated by auxiliary laser 
fields. 

In the next section, we study the properties of the op- 
tical memory for a finite optical depth when the atomic 
ensemble reemits the light pulse in forward and back- 
ward directions. We also look at the impact of the ini- 
tial distribution bandwidth and wc establish the optimal 
broadening of the initial distribution with respect to the 
memory efficiency for a given pulse. 



III. GENERAL SOLUTION 



We solve the equations of motion by generalizing the 
technique presented in ref. [2^ . To clearly distin- 
guish the contribution of the initial distribution of atoms 
Go(Ao) and of the broadened distribution G'(A') associ- 
ated to a single absorption line, we introduce the atomic 
operator 



CTij(z,i;A',Ao) 



1 



7V(A' + Ao,z: 



Af(A' + Ao,2) 
n=l 



(9) 

where A^(A' -I- Ao,z) is the number of atoms with ini- 
tial detuning around Ao broadened to A' + Aq and with 
position around z. When the inhomogeneous broaden- 
ing is reversed, the detuning is changed from A' + Aq 
to — A' -f Ao such that the atomic properties are de- 
scribed by the operator (Tij(z, t\ — A', Aq). As in the pre- 
vious section, the atomic operator is decomposed into two 
counter-propagating contributions. We are interested in 
the evolution of an incoming light pulse propagating in 
the forward direction through the atomic ensemble 



( a a \ 

Kdt^'Tz) 

/+00 
dAodA'G(Ao)G'(A')af (z, t; A', Aq), 
-00 



(10a) 
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— af(z,t;A',Ao) 



(10b) 



-i{Ao + A')(Tf (z, t; A', Ao) + tpEi''{z, t), 

Initially, there is no atomic excitation such that (J{{z,t 
— oo; A', Aq) — 0. Taking into account this initial condi- 
tion, Eqn. (jlObp can be solved as 



CTf(z,i;A',Ao) 



(11) 

We are interested in the regime for which the storage 
duration T is longer than the pulse duration r. We thus 
consider that at the position z = 0, the incoming light 
pulse E™{z — 0) is centered around — T/2. The inhomo- 
geneous broadening is reversed at time t = Q after a time 
sufficiently long with respect to the pulse duration r. The 
Fourier transform of the incoming light pulse is defined 

by 



(12) 



Ef{z,uj) := / dt e"^'El'^{z,t). 



The upper bound is equal to since the incoming light 
pulse is not defined anymore at positive times. Further- 
more, wc consider the regime in which the broadened 
distribution bandwidth 7 is larger than the inverse of the 
storage duration 1 /T. These considerations allow us to in- 
troduce the Fourier transform of the incoming light pulse 
in Eqn. pOap in which Cf is replaced by its expression 

m 



— - — +flH{u:)]Ef\z,uj)=0. 



(13) 



H{lo) is defined by 



H{uj) j dxie"^"" x 



dAodA'Go(Ao)G'(A')e-^('^'+'^«)") 



and 77 := g^Np^ /c. The solution of Eqn. ([T^ is given by 



Er{z,Lo)^E^[0,uj)e 



iujz / c —r]H{uj)z 



(14) 



at any position inside the medium ^ z ^ L. The solu- 
tions pT|) and describes the absorption of the incom- 
ing light pulse by the atomic ensemble. The absorption 
coefficient a{Lo) defined from the relation |£^™(2;,a;)p = 
e""('^)^|£;^"(0,cj)|2 is given by 



a{Lo) 



2g2iVp2Re(H(c^)) 



(15) 



It is thus proportional to the coupling between the light 
pulse and the atomic ensemble and depends on the 
atomic distribution via the function H{uj). 



In what follows, we give the expressions of the outgoing 
light pulse reemitted by the atomic ensemble by distin- 
guishing the pulses reemitted in the controlled backward 
and forward directions. 



A. Complete protocol: Emission in backward 
direction 



The light pulse is partially absorbed such that a part 
of the light excitation is transfered into the atomic en- 
semble. At time t = 0, the inhomogeneous broadening 
is reversed meaning that the detuning Aq -I- A' changes 
to Aq — A'. Furthermore, the phase shift e^'^^z/c jg g^p_ 
plied to the atomic ensemble such that the system evolves 
backward in space. We assume that at t = the non- 
absorbed part of the light has left the medium so that 
the field is zero, i.e. the excitation is purely atomic. The 
system is described by the following set of equations 



(16a) 



i(3 



/+00 
dAodA'Go(Ao)G(-A')ab(z, t; -A', Aq), 
-00 



d_ 

dt 



crb(z,i;-A',Ao) = 



(16b) 



-z(Ao - A')ab(z, t; -A', Ao) + ipE^^^z, t) 
with the initial condition 

ii;r'(^, i = o) = o, 

/•o 

crb(z, t = 0; -A', Ao) =ip ds e'^^ +'^°>EY'{z, s). 

<J —00 

The solution of Eqn. (jl6b[l is given by 

ab(z, <; -A', Ao) = a,,{z, 0; -A', Ao) + ds[e^('^'-^«)(*- 
X {ip E°^\z, s) + liA' - Ao) ab(z, 0; -A', Ao)) ] . 



Plugging the previous expression into Eqn. p6ap . we get 
the following propagation equation 



§-^ + '^^vFH]Er(z,u;) 



(17) 



r, / dAoGo(Ao)J(w;Ao)£;^"(z,-a; + 2Ao) 

J —oo 

for the Fourier transform of the outgoing pulse 

r+oc 

E^''\z,uj) -.^ dt e''^'E^''\z,t). (18) 

Jo 

The lower bound is equal to since the outgoing light 
pulse, centered around T/2, is not defined at negative 
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times. We have introduced the functions 



forward direction is given by 



J(w;Ao) 
and 



dx 



dA'G'(-A')e'('^'-'^«)^e 



+ 00 



dA 



odA'Go(Ao)G"(-A')e'(^'-^'''") . 



The solution of Eqn. (fT7|) estabhshes the connection be- 
tween the output and the input hght pulses when they 
propagate in opposed directions 

r + ca 

K'^'iz, Lj)^-v dAoGo(Ao)e-^-^/^ x 



^2iAoL/c^-r,LH{-u+2Ao)^v{z-L)F{Lo) 



g2iAoz/Cg-77xif(-w+2Ao) 



J(w;Ao) 



2iAo/c - riH{-uj + 2Ao) - ■qF{ijj) 



^j"(0,-c^ + 2Ao) 



(19) 



In the next subsection, we study the simplified protocol in 
which the position-dependent phase is not applied to the 
atomic ensemble. We find the expression of the output 
pulse when it propagates in the forward direction. 



B. Protocol without phase shift: Emission in 
forward direction 



If at time t = 0, the additional phase e^™"^/'^ is not 
applied to the atomic ensemble, the system evolves for- 
ward in space and its evolution is given by the set of 
equations (|10p in which A' has to be changed to —A'. 
As above, the light pulse is partially absorbed and the 
absorbed light excitation is transfered to the atoms. The 
initial condition is thus 

£;f™*(z,t = 0) = 0, 

ai{z, t = 0; -A', Aq) = ip / ds e'^^ +^°^'E^{z, s). 



The equations of motion are resolved as previously. The 
expression of the outgoing light pulse propagating in the 



+00 



E'i^'{z,uj) = -77Z / dAoGo(Ao)J(w; Ao) x 

<J —OQ 

. , / F{u;) -H{^iu + 2 Ao) . u; - Ao\ 
smhc rjz iz | x 



Ao F{uj) + H{-oj + 2Aoy 
exp [ iz rjz I X 

^^"(0,-c^ + 2Ao). 



(20) 



The sinhc denotes the hyperbolic sinus cardinal function 
(sinhc(a;) — sinh(x)/x). 

The equations and (|^ establish the expression of 
the light pulse reemitted by the atomic ensemble in both 
backward and forward directions as a function of the in- 
put light pulse characteristics and of the properties of 
the atomic medium. They are valid for a general initial 
distribution broadened to an arbitrary distribution. In 
the next sections, they are analyzed in order to deduce 
some characteristics of the memory for specific atomic 
distributions. 



IV. FINITE OPTICAL DEPTH 

In the framework of a practical realization, it is 
necessary to know the properties of the memory for a 
finite optical depth. This topic has also been addressed 
for the CRIB protocol in Refs. pi, 12^1 and for the 
slow-light based memories in Refs. 25, 26]. 

The analysis is done using a simple model in which 
the initial distribution is reduced to a single narrow 
absorption line Go(Ao) = <^(Ao). The initial absorp- 
tion line is broadened to be constant on the interval 
[—7/2,7/2] and null elsewhere such that it satisfies the 
condition /+^dA'G'(A') ^ dAG(A) = 1. In the 
regime where the spectral pulse bandwidth is smaller 
than the inhomogeneous broadening F <C 7, we have 
H{uj) w F{uj) « J{u!; 0)/2 w tt/j and the outgoing pulse 
propagating in the backward direction (fTO)) takes the fol- 
lowing form 



£;r'(o,t) = -(1 -e-"^)£;^"(o,-<). 



(21) 



at the position z — 0. The absorption coefficient is given 
by a = 2TTg^Np'^/"fc. The memory efficiency, defined by 



Eff := 



/dt^lE™*(t^)l2 

/dt^lEi"(cj)|2 



(22) 



corresponds to the probability to retrieve an absorbed 
photon. The memory efficiency depends on the optical 
depth aL and as can be seen in Fig. O it reaches 100% for 
large enough optical depth. This limit of a large optical 
depth corresponds to the result presented in Ref . [ig] . 
Under similar conditions, when no phase shift is ap- 
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FIG. 2: (Color online) Efficiency of the light pulse reemission 
in forward (full blue line) and backward (dashed black line) 
direction as a function of the optical depth for a constant 
broadening. The efficiencies vary as (aL)^ e~°'^ for the for- 
ward protocol and as (1 — e~"^)^ for the backward protocol. 



where /9o(Ao) = 7VGo(Ao) is the initial density of atoms. 
uqL corresponds to the initial optical depth in the ab- 
sence of broadening. For a given medium , the optical 
depth aL varies with the ratio of the initial distribution 
bandwidth on the broadened distribution bandwidth. In 
a first subsection, we show that the initial distribution 
bandwidth also limits the storage duration. We thus 
prove the existence of a trade-ofF between the memory 
efficiency and the storage duration. In a second subsec- 
tion, we look at the optimal broadening for a given pulse 
and for a fixed initial distribution. 



plied, the expression of the output light pulse propagat- 
ing in the forward direction (j20p is given at the position 
z = L by 



£:f°"*(L,tj) = -aLe 



~aL/2 sin(^i^/c) ~ 



bjL/c 



ETiO.-uj). (23) 



The sinus cardinal function induces a distortion of the 
light for large pulse bandwidth F with respect to c/L. 
For a pulse bandwidth smaller than c/L, the expression 
for the forward pulse is reduced to 



E^''\L,t) = -aLe 



-aL/2 



i?r(o,-i). 



(24) 



In this regime, the interaction with the atomic ensemble 
does not induce distortion neither in the forward direc- 
tion nor in the backward direction. For small optical 
depth, the memory efficiency varies quadratically as a 
function of the optical depth for both forward and back- 
ward protocols. However, in the forward direction, the 
pulse is reabsorbed by the atoms when the optical depth 
increases and the memory efficiency reaches a maximum 
of 54% at the optical depth aL — 2. Here we have sup- 
posed that the optical depth aL is constant for the whole 
spectral width of the pulse. If aL is lower for certain 
spectral components of the pulses, the memory efficiency 
is lower. 



V. ATOMIC DISTRIBUTION 

In the previous section, we have studied the depen- 
dence of the memory efficiency on the optical depth. 
For an initial distribution with characteristic bandwidth 
Go(Ao) ~ 7o smaller than the spectral pulse bandwidth 
F and broadened to a distribution bandwidth G(A) « 7, 
the optical depth can be written as 



(25) 



A. Memory efficiency versus storage duration 



To analyze the role of the initial distribution band- 
width, we consider the simple case in which the broad- 
ened distribution is constant in a given interval and null 
elsewhere 



1 



G'(A') = -9{A' 



1 



l)eC- 
2' 4 



A') 



(26) 



{d being the Heaviside function) such that in the regime 
F ^ 7 and 70 <^ 7, we have J(a;,Ao)/2 w H{i^) ~ 
F{uj) « n/'-f. We discuss the results obtained from the 
backward process but the conclusions are applicable for 
the forward process as well. Under the condition 70 ^ 
c/L, the output light pulse emitted in the backward di- 
rection (|19p is given by the following expression 



E°''\o, t) = -(1 - e-"-^)£;;"(o, -t) X 

+ 00 

dAoGo(Ao)e-2'^«* 



(27) 



at the position z = 0. Since the output pulse is cen- 
tered around T/2, we are interested in times t around 
r/2. Comparing Eqns. (|2T|) and (|27p. it clearly appears 
that the output pulse is multiplied by the Fourier trans- 
form of the initial distribution. For an initial distribu- 
tion bandwidth 70, the storage duration T is thus limited 
by 1/70- This has already been observed experimentally 
as reported in ref. [17]. The bandwidth of the initial 
distribution thus constitutes a limitation for the storage 
duration T70 <^ 1. Since the initial distribution width 
also affects the optical depth (see Eqn. (|25|) ) and thus 
the memory efficiency, there exits a trade-off between 
a long storage duration and a large memory efficiency. 
Note that this is true for a two-level system. If other 
long-living levels are available, the storage time can be 
increased significantly by transferring the excitation to 
a second long-living state, e.g. a ground state hyperfine 
level [13. 
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B. Optimal broadening 

We now consider that the initial distribution has been 
chosen. We define the effective width of the initial dis- 
tribution 



glNp^L g2pVo(0)7oi 



(28) 



which depends on the medium properties together with 
the initial distribution and thus has a fixed value. The 
optical depth after broadening is given by 



(29) 



For a given pulse shape, we look for the optimal width of 
the broadened distribution with respect to the memory 
efficiency. 

We first fix the shape of the atomic distributions to be 
the same as the Lorentzian shape of the light pulse 



Go(Ao) - 
G'(A') = 



1 



27r r2/4- 
70 1 



^0 



7 1 

2iT 72/4 + A'2' 



(30) 

(31) 
(32) 



For a given value of the effective width of the initial dis- 
tribution v, we calculate the memory efficiency by vary- 
ing the broadened distribution width with respect to the 
pulse bandwidth. In the limit 70 — > 0, po(0) +00 such 
that V has a finite value, the functions J(w; 0), F{uj) and 
Hiuo) take the simple forms 



J(c^;0) 



7 



72/4-1-^^2' 



The output pulse HI])- ([201) is given by 



(33) 



^r(o,c.) = - 



27r(r2/4-Fcj2) 



1 — exp — 



72/4 + Cj2 

(34) 



27r(r2/4-htj2) ^2/4 



exp 



1/7 



2(72/4 -ftj2) 



smc 



^2 -f -^2 



(35) 



when propagating in the backward and forward (with the 
hypothesis F smaller than c/L) directions respectively. 
In Fig. [3] (top), we take v = 2T and we plot the efficiency 
as a function of the broadened distribution bandwidth 7. 
For large 7, the optical depth ((25)) which is inversely pro- 
portional to 7, is small. As shown in Fig. [3] (top), the 
memory efficiency for the backward protocol is close to 
the efficiency of the forward protocol. For smaller val- 




FIG. 3: (Color online) Efficiency of the light pulse reemission 
for forward (full blue line) and backward (dashed black line) 
protocol as a function of the broadened distribution band- 
width 7 (in units of the pulse bandwidth F) for effective 
widths of the initial distribution u = 2T (top) and v — O.OSF 
(bottom). 



ues of 7, the efficiency increases and reaches a maximum 
for the forward process caused by the reabsorption while 
the efficiency of the backward protocol continues to grow 
until a value close to 100%. When 7 tends to zero, the 
bandwidth of the broadened distribution is too small to 
absorb the light pulse and the memory efficiency tends to 
zero for both backward and forward processes. In both 
cases, however, the optimum efficiency occurs when the 
broadened distribution bandwidth and the light band- 
width are of the same order. 

Taking a smaller value for v — 0.05F, we plot in Fig. [3] 
(bottom) the emission efficiency versus the distribution 
bandwidth. The efficiency curves exhibit similar shapes 
as in the large v case except that their maxima reach 
smaller values. These maxima are obtained for a nar- 
rower broadened distribution bandwidth than the pulse 
bandwidth. 

One can understand these results in the following way: 
To insure a high efficiency of the CRIB protocol, the op- 
tical depth after broadening has to satisfy aL > 1. From 
Eqn. ([^5]) , we deduce that the effective width of the ini- 
tial distribution has to be at least of the order of the 
broadened distribution bandwidth: v ^ "f- Thus, for a 
pulse bandwidth F w z^, one can broaden the initial dis- 
tribution until 7 « F and thus achieve high efficiency. 
If F ^ z/, it is more advantageous to choose a smaller 
broadening than the pulse width 7 « z^ <C F. In this 
case, however, one must expect strong distortion of the 
output pulse due to filtering effects. 
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C. Broadening shape 



pulse generated in the forward direction (120]) is of the 
form 




TTlOUt 



(42) 



FIG. 4: (Color online) Comparison of memory efficiencies for 
a spectral distribution of atoms and light pulse with identical 
shapes (forward: full blue line, backward: dashed black line) 
and with different shapes (forward: dashed-dotted red line, 
backward: dotted red line) as a function of the inhomogeneous 
bandwidth (in units of the pulse bandwidth). The effective 
width of the initial distribution v is taken equal to 2r. 

Finally, we investigate how the shape of the broadened 
distribution influences the memory efficiency. We keep 
the Lorentzian shape of the light pulse but we take a 
constant distribution in a given interval 



-Er(0,w) 



r 



1 



2tt r2/4 + w2' 



(36) 



Go(Ao) = —9{Ao + 7o/2)0o(7o/2 - Ao), (37) 
70 

G'(A') = ie(A' + 7/2)0(7/2 -A'). (38) 

7 

For a fixed value of the effective width of the initial dis- 
tribution v, we calculate the memory efficiency and we 
compare it to the efficiency obtained for the Lorentzian 
shapes of Eqns. ([30ll32l) . 

In the limit 70 0, the functions J, F and H take the 
following forms 



J{iu;0) 



2tt 



(a; + 7/2)0(7/2-^), 



(39) 



FH-i^(-r = ^ + -iog^, (40) 

2 7 7 — 2u; 

such that the light pulse generated in the backward di- 
rection (|19p is given by 



r 



1 - 



1 



(41) 



27r V " JV^/A + uj"^ 
for Ld e [—7/2,7/2] and elsewhere. For F <C c/L. the 



7 r2/4 + cj2 

for to € [—7/2,7/2]. From these expressions, we plot (see 
Fig. 2]) the memory efficiency. We compare it to the 
efficiency we found for identical broadened distribution 
and light pulse shapes. Globally, the relative shapes of 
the broadened distribution and of the light pulse does not 
modify the memory efhciency. Locally, for small distribu- 
tion bandwidths, the memory efficiency is a bit larger for 
identical shapes. However, when the distribution band- 
width increases, the efficiency is higher when these shapes 
are different. 



VI. CONCLUSION 

We have obtained the explicit solution of the equations 
of motion for the system atoms plus light in the weak 
excitation regime, making it possible to gain insight 
into the dependence of the memory efficiency on the 
optical depth, and on the width and shape of the atomic 
spectral distributions. Furthermore, we introduced a 
simplified CRIB-based memory protocol, which does 
not require the use of any additional laser field, and 
showed that its efficiency is limited to 54%. For both 
the complete and the simplified CRIB protocols, the 
interaction with the atoms does not induce distortion of 
the stored light pulse. We have shown that the storage 
duration is limited by the initial spectral distribution of 
the atoms. This distribution also defines the absorption 
of the medium which determines the memory efficiency. 
There thus exists a trade-off between the storage du- 
ration and the memory efhciency. For a given pulse 
shape, the optimal broadening of the initial distribution 
is of the order of the effective bandwidth of the initial 
atomic distribution, defined as the product of the atomic 
density and the width of the initial atomic distribution. 
We have also shown that the shape of the broadened 
distribution does not dramatically change the memory 
efficiency. 
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